Large-scale Monte Carlo simulations of the isotropic three-dimensional Heisenberg 

spin glass 
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We study the Heisenberg spin glass by large-scale Monte Carlo simulations for sizes up to 32 3 , 
down to temperatures below the transition temperature claimed in earlier work. The data for the 
larger sizes show more marginal behavior than that for the smaller sizes, indicating the lower critical 
dimension is close to, and possibly equal to three. We find that the spins and chiralities behave in 
a quite similar manner. 
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I. INTRODUCTION 

Following the convincing numerical work of Ballesteros 
et alii there has been little doubt that Ising spin glasses 
in three dimensions have a finite temperature transition. 
In this paper we shall study a related model for which 
the existence of a finite temperature transition is still con- 
troversial: the isotropic Heisenberg spin glass, which is 
composed of classical spins with three components. Early 
wor k 2 ' 3 ' indicated a zero temperature transition, or pos- 
sibly a transition at a very low but non-zero temperature. 
However, following the pioneering work of Villain 5 , which 
emphasized the role of "chiralities" (Ising-like variables 
which describe the handedness of the non-collinear spin 
structures), Kawamura and Tanemura 6 proposed, for the 
XY case (which have two-component spins) , that the spin 
glass transition only occurs at Tsg = and that a chi- 
ral glass transition occurs at a finite temperature Tcg- 
This scenario requires that spins and chiralities decou- 
ple at long length scales. Kawamura and collaborators 
subsequently proposed2>2»£ that this "spin-chirality de- 
coupling" scenario also holds for Heisenberg spin glasses. 

However, the absence of a spin glass transition in 
Heisenberg spin glasses has been challenged by Mat- 
subara et ah 10 ' 11 , and Nakamura and Endohi 2 - who ar- 
gued that the spins and chiralities order at the same 
low but finite temperature. Recently Picco and Ritort 13 
also claimed a finite Tsg and inferred that probably 
Tsg — Tcg, though they did not investigate the chi- 
ralities directly. In earlier workii, referred to as LY, we 
studied spin and chiral correlations on an equal footing, 
using the method of analysis that was the most success- 
ful for the Ising spin glassi, namely finite-size scaling of 
the correlation length. Considering a modest range of 
sizes, L < 12, we found that the behavior of spins and 
chiralities was quite similar and they both had a finite 
temperature transition, apparently at the same temper- 
ature. 

However, quite recently Campos et alii 5 - were able to 



study larger sizes than LY, up to L — 32. They agreed 
with LY that there is a single transition at which both 
spins and chiralities order, but they also found evidence 
for crossover, at the largest sizes, to a "marginal" be- 
havior, reminiscent of that at the Kosterlitz-Thouless- 
Berezenski i 16 ' 17 (KTB) transition in the two-dimensional 
XY ferromagnet where there is a finite transition tem- 
perature T c but no long-range order for T < T c . In 
fact the region below T c is a line of critical points in 
the KTB theory. If a line of critical points also exists in 
the three-dimensional Heisenberg spin glass, then d = 3 
is the lower critical dimension di, below which there is 
no transition. However, given numerical uncertainties, 
Campos et al. cannot rule out the possibility that di is 
slightly less than three, in which case there is spin glass 
order below Tsg- 

Hukushima and Kawamura 9 studied sizes up to L = 
20 and found more marginal behavior when comparing 
L = 16 and 20, than for the smaller sizes. However, they 
argued that this effect is greater for the spin correlation 
length than for the chiral correlation length, and hence 
concluded that, while Tcg is finite, Tsg is zero or possi- 
bly non-zero but less than Tgg-, i-e. there is spin-chirality 
decoupling. 

In this paper we perform Monte Carlo simulations of 
the Heisenberg spin glass, along the lines of LY, but for 
larger sizes. Our main motivation is to investigate the 
claim of Campos et alii 5 - that there is a line of critical 
points for T < Tsg- In order to test whether there is a 
critical line, as proposed by Campos et al, or the usual 
situation of a single critical point at Tsg, it is necessary 
to investigate the behavior of the system below the esti- 
mated Tsg- Campos et al. were not able to do this and 
the evidence for the critical line was based on estimat- 
ing corrections to scaling at Tsg- This is rather indirect, 
and perhaps not very reliable because the range of sizes 
and quality of the data, are not sufficient to disentangle 
the various corrections to scaling unambiguously. Very 
recently, the analysis of Campos et al. has been criti- 
cized by Campbell and Kawamura^. Here, we do not 
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rely on corrections to scaling but obtain data below T$g 
and so directly find evidence for more marginal behavior 
at larger sizes. 

As we shall see, the transition temperature in the 
Heisenberg spin glass is very low and so it is surpris- 
ing to us that Campos et al^ did not use the technique 
of "parallel tempering" 19 ! 20 , which is the commonly used 
approach to speed up simulations of spin glasses, espe- 
cially at very low temperatures. Equilibrating lattices as 
large as 32 3 is very challenging, especially in the absence 
of parallel tempering. A second motivation of our study 
is therefore to use parallel tempering and to combine this 
with a very useful equilibration test described in the next 
section, to ensure that the data are equilibrated. 

From our results, we see a strong crossover to much 
more marginal behavior for sizes L > 24, in agreement 
with Campos et al.— . Whether there is a KTB-like crit- 
ical line as proposed by Campos et al., is, however, un- 
clear. The lower critical dimension, di seems to be close 
to, or possibly equal to, three. The behavior of the spins 
and chiralities is rather similar, and so, in contrast to 
Kawamura and collaborators 9 *^, we do not feel that our 
data supports spin-chirality decoupling. 

II. MODEL AND ANALYSIS 

We use the standard Edwards-Anderson spin glass 
model 

^ = — JijSi ■ Sj, (1) 

where the are 3-component classical vectors of unit 
length at the sites of a simple cubic lattice, and the Jy 
are nearest neighbor interactions with a Gaussian dis- 
tribution with zero mean and standard deviation unity. 
Periodic boundary conditions are applied on lattices with 
N = L 3 spins. 

The spin glass order parameter, g A " / (k), at wave vector 
k, is defined to be 

^(k) = l^^ 1 ^r (2) e lk ' R % (2) 

i 

where /i and v are spin components, and "(1)" and "(2)" 
denote two identical copies of the system with the same 
interactions. From this we determine the wave vector 
dependent spin glass susceptibility xsgO*-) by 

X5G(k) = A^[(|-r(k)| 2 )] av , (3) 

where (• • •) denotes a thermal average and [• • -] av de- 
notes an average over disorder. The spin glass correlation 
length is then determinedi^, from 

1 ( Xsg(0) V /2 ; (4) 

2sin(fc min /2) \XSG(kmin) / 



where k min = (2tt/L)(1, 0, 0). 

For the Heisenberg spin glass, Kawamura 7 defines the 
local chirality in terms of three spins on a line as follows: 

K i = Si-^ ■ 8^ x Si— p,. (5) 

The chiral glass susceptibility is then given by 

X £ G (k)=A[(|^(k)| 2 )] av , (6) 

where the chiral overlap (k) is given by 

^( k ) = ^E< (1) < (2)e ' ik ' Rl - ( 7 ) 

i 

We define the chiral correlation lengths ^ L by 

^ L 2sin(fc min /2) Uc G ( k min) J ' 

in which xcgO*- — 0) is independent of fi. Note that 
£c l wn l' m general, be different for p, along k min (the 
x direction) and perpendicular to k, though we expect 
that this difference will vanish for large sizes. We denote 
these two lengths by Q L and respectively. 

To equilibrate the system in as small a number of 
sweeps as possible, with the minimum amount of CPU 
time we perform three types of Monte Carlo move: 

1. "Microcanonical" sweeps^, 

(also known as "over-relaxation" sweeps). We 
sweep sequentially through the lattice, and, at each 
site, compute the local field on the spin, = 
^2 j Jij . The new value for the spin on site i 
is taken to be its old value reflected about H, i.e. 

S; : = -S, + 2 S ^H l , (9) 

see Fig, [TJ These sweeps are microcanonical be- 
cause they preserve energy. They are very fast 
because the operations are simple and no random 
numbers are needed. For reasons that are not fully 
understood, it also seems that they "stir up" the 
spin configuration very efficiently 15 and the system 
equilibrates faster than if one only uses "hcatbath" 
updates, described next. 

2. "Heatbath" sweeps! 

Since the microcanonical sweeps conserve energy 
they do not equilibrate the system. We therefore 
also include some heatbath sweeps since these do 
change the energy, typically doing one after every 
10 microcanonical sweeps. As for the microcanoni- 
cal case, we sweep sequentially through the lattice. 
Referring to Fig.[TJ we take the direction of H,;, the 
local field on site, i to be the polar axis, denote the 
polar angle as 9, and define the azimuthal angle <p 
such that (p — for the old spin direction. The new 
spin direction is characterized by angles 9 and 
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FIG. 1: Hi is the local field on site i due to its neighbors. The 
spin at i is initially in direction Si. In a microcanonical (over- 
relaxation) move, the spin is reflected about Hi according to 
Eq. ©, and so ends up in direction 



</>, relative to H^, as follows. The energy does not 
depend on the azimuthal angle, and so <fi = 27rri, 
where n is a random number chosen uniformly be- 
tween and 1. The polar angle is chosen such that, 
after the move, the spin is in local equilibrium with 
respect to the local field, i.e. if x = cos 8, then 



P(x) 



2 sinh 0Hi 



(10) 



where (3 = l/T. To determine x with this proba- 
bility, the procedure^ is to equate the cumulative 
distribution 



Q(x) 



P(x') dx 



-0Hi 



(11) 



to a second random number r 2 (also with a uniform 
distribution between and 1) so 



1 



ln[l + r 2 (e^'-l)] -1. 



(12) 



It is then necessary to convert the new spin direc- 
tion back to cartesian coordinates. Denoting 
the polar and azimuthal angles of H by 9h and 4>h 
relative the cartesian reference frame, and remem- 
bering that and cj> are relative to H, we have 2 ^ 



S' x = S" cos 4>h — Sy sin 4> H 
= S" sin <j> H + S" cos 4>h 



S' z 
where 

S ,. 

s" 



cos 6 cos Oh — sm 6 sm Oh cos < 



cos sm Uh + sm cos Oh cos ( 



(13) 
(14) 
(15) 



(16) 
(17) 



We see that the calculations in the heatbath moves 
are quite involved, which is why we do mainly mi- 
crocanonical moves, just including some heatbath 
moves to change the energy and thereby ensure the 



algorithm is ergodic. Note, though, that the accep- 
tance probability for the heat bath moves, and also 
for the microcanonical moves, is unity, so no moves 
are wasted. 

"Parallel tempering" sweeps. 

At low temperatures spin glasses are easily trapped 
in minima (valleys) of the free energy. In order 
to ensure that the system visits different minima 
with the correct Boltzmann weight during the time 
of the simulation we use the method of parallel 
temperingAS»22'. One takes Nt copies of the system 
with the same bonds but at a range of different tem- 
peratures. The minimum temperature, T m i n = T\, 
is the low temperature where one wants to investi- 
gate the system (below T$g in °ur case), and the 
maximum, T max = 7V T , is high enough that the 
the system equilibrates very fast (well above T$g 
in our case). A parallel tempering sweep consists 
of swapping the temperatures of the spin config- 
urations at a pair of neighboring temperatures, Ti 
and Ti+i, for i = 1, 2, • • • , Tjy T _i with a probability 
that satisfies the detailed balance condition. The 
Metropolis probability for this is^ 



P(T swap) 



exp(-A^A£:), (if A/3 AS > 0), 
1, (otherwise), 

(18) 

where A/3 = 1/Ti-1/T i+1 and AE = Ei~E i+1 , in 
which Ei is the energy of the copy at temperature 
Ti. In this way, a given set of spins (i.e. a copy) per- 
forms a random walk in temperature space. Sup- 
pose that at some time in the simulation a copy is 
trapped in a valley at low-T. Later on it will reach 
a high temperature where it randomizes quickly, so 
that when, still later, it is again at a low temper- 
ature, there is no reason for it to be in the same 
valley that it was in before. We do one sweep of 
temperature swaps after every ten microcanonical 
sweeps. 

Table U gives the parameters of the simulations. It 
will be seen that the number of temperatures is very 
large, and it is instructive to discuss the reason for 
this. The difference between two neighboring tem- 
peratures, AT, must be sufficiently small that there 
is an overlap between the energy distributions at 
those temperatures, so there is a reasonable prob- 
ability that the same spin configuration occurs for 
both temperatures. Otherwise, the probability for 
accepting the temperature swap, Eq. (fl"8|) , will be 
very small. Relating the width of the temperature 
distribution to the specific heat in the normal way, 
one finds that 



5T 



< 



1 



Vcn' 



(19) 



where C is the specific heat per spin. For the 
Heisenberg spin glass C tends to a constant at low- 
T because of the Gaussian (spinwave) fluctuations 
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TABLE I: Parameters of the simulations. N samp is the num- 
ber of samples, -/V equ ii is the number of microcanonical Monte 
Carlo sweeps for equilibration for each of the 2Nt replicas for 
a single sample, and iV meas is the number of microcanonical 
sweeps for measurement. The number of heatbath sweeps is 
equal to 10% of the number of microcanonical sweeps. Tmm 
and T max are the lowest and highest temperatures simulated, 
and Nt is the number of temperatures used in the parallel 
tempering. 
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N 

1 ' samp 


Acquil 


1 v mcas 


T ■ 

1 mm 


T 

1 max 


N T 


4 


500 


5 x 10 3 


10 4 


0.0400 


0.96 


40 


6 


500 


6 x 10 4 


1.2 x 10 5 


0.0400 


0.96 


40 


8 


536 


8 x 10 4 


1.6 x 10 5 


0.0400 


0.96 


57 


12 


204 


8 x 10 5 


1.6 x 10 6 


0.0400 


0.61 


88 


16 


202 


8 x 10 5 


1.6 x 10 6 


0.1015 


0.49 


73 


24 


160 


3.2 x 10 5 


6.4 x 10 5 


0.1200 


0.49 


122 


32 


56 


9.6 x 10 5 


1.28 x 10 6 


0.1210 


0.40 


120 



about the local equilibrium positions. Hence, for a 
given size, we choose temperatures which decrease 
in a geometric manner. The required number of 
temperatures to cover a fixed range of T increases 
proportional to y/N, and so, for large N, many 
copies are needed to cover even a factor of 2 in 
T. However, for the Ising spin glass with Gaus- 
sian interactions, the specific heat tends to zero as 
T — > (presumably linearly in T). Hence much 
bigger steps in T can be taken at low-T than for 
the Heisenberg case, leading to fewer temperatures 
being needed. 

To test for equilibration^ we require that data satisfy 
the relation^ 

qi-q s = - z TU, (20) 

which is valid for a Gaussian bond distribution. Here 
U = -E(i.j) Jij{Si-Sj)]av is the average energy per spin, 
m = (1M) E(, J ) [(Si • S,) 2 ] av is the "link overlap", q s = 
(Wb) E(ij>K(Si ' S 3 ) 2 )]av, N b = (z/2)N is the number 
of nearest neighbor bonds, and z (= 6 here) is the lattice 
coordination number. Equation (|20[) is easily derived by 
integrating by parts the expression for the average energy 
with respect to Jjj, noting that the average [• • -] av is over 
a Gaussian function of the J^'s. 

The spins are initialized in random directions so the 
energy, the RHS of Eq. (|2"D|) , is initially close to zero and 
decreases, presumably monotonically, to its equilibrium 
value as the length of the simulation increases. Hence 
the RHS of Eq. ([2H)l will be too large if the simulation is 
too short to equilibrate the system. On the other hand, 
the LHS of Eq. (f2"0"|) , qi will be too small if the simula- 
tion is too short because it starts off close to zero and 
then increases with MC time as the two replicas start 
to find the same local minima. The quantity q s will be 
less dependent on Monte Carlo time than qi since it is a 
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FIG. 2: (Color online) Equilibration plot, testing Eq. ([20|) for 
L = 24 at the T = 0.144. It is seen that the data for 2TU/z 
come together at about 3 x 10 5 total sweeps (equilibration 
plus measurement) and then stay at their common value in- 
dicating that equilibration has been achieved. The lines are 
guides to the eye. It is seen that the energy comes close to 
its equilibrium value very quickly, whereas qi — q 3 takes much 
longer. 



local variable for a single replica. (For the Ising case it is 
just a constant.) Hence if the simulation is too short the 
LHS of Eq. (|2H|) will be too low. In other words, the two 
sides of Eq. (|20p approach the common equilibrium value 
from opposite directions as the length of the simulation 
increases. Only if Eq. (|20p is satisfied within small error 
bars do we accept the results of a simulation. 

Figure shows a test to verify that Eq. PD)) is sat- 
isfied at long times. For the parameters used, L = 
24, T = 0.144, this occurs when the total number of 
sweeps (N sw = N cqui i + N mcas ) is about 3 x 10 5 . Figures 
[3] and H] show that the spin and chiral correlation lengths 
appear to become independent of N sw , and hence are pre- 
sumably equilibrated, when N sw is larger than this same 
value. Hence, it appears that when Eq. (|20p is satisfied 
to high precision, the data for the correlation lengths is 
equilibrated. 

With the number of sweeps shown in Table QJ Eq. (|2"0|) 
was satisfied for all sizes and temperatures. The error 
bars are made sufficiently small by averaging over a large 
number of samples. We are simulating system sizes which 
are very large by spin glass standards (up to N = 32 3 ), 
so it is crucial to have a stringent test like this for equi- 
libration. 

Since £l / L is dimensionless it has the finite size scaling 
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FIG. 3: (Color online) A plot of £l/L as a function of the total 
number of sweeps for L — 24 at T = 0.144. It is seen the data 
flattens off at around 3 x 10 5 sweeps, the value where the 
two sets of data in Fig. [2] start to agree. This indicates that 
when the data in Fig. [2] agree within high precision, i.e. when 
Eq. (|20|l is satisfied, the correlation length has reached its 
equilibrium value. The line is a guide to the eye. 



L 



X 



(lV"(T- 



Tsg) 



(21) 



where v is the correlation length exponent. Note that 
there is no power of L multiplying the scaling function 
X, as there would be for a quantity with dimensions. 
There are analogous expressions for the chiral correla- 
tion lengths. From Eq. (|2ip it follows that the data for 
£l/L for different sizes come together at T = Tsg- In 
addition, they are also expected to splay out again on the 
low-T side if there is spin glass order below Tsg- In a 
marginal situation, with a line of critical points as in the 
KTB transition, the data for different sizes would come 
together at Tsg and then stick together at lower T, see, 
for example, Fig. 3 of Ref. QJ 



III. RESULTS 



FIG. 4: (Color online) A plot of Q L /L as a function of the 
total number of sweeps for L = 24 at T = 0.144. As for Fig. [3] 
it is seen the data flattens off at around 3 x 10 5 sweeps, the 
value where the two sets of data in Fig. [2] start to agree. The 
line is a guide to the eye. 
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We studied sizes from L — 4 to L = 32, as shown in 
Table HI The CPU time involved to get this data is about 
15 Mac G5 CPU years. 

We start with the data for shown in Fig. [51 As 

was found earlier by LY, the smaller sizes show a clear 
intersection, and also splay out at lower temperatures 



FIG. 5: (Color online) Data for £l/L, the spin glass corre- 
lation length divided by system size, as a function of T for 
different system sizes. Note that there are very many data 
points for the larger sizes. This is because a large number of 
temperatures are needed for the parallel tempering algorithm, 
as discussed in Sec. [TTJ 
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FIG. 6: (Color online) Enlarged view of a region of Fig. [5] 



FIG. 8: (Color online) Enlarged view of a region of Fig. [7] 
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FIG. 7: (Color online) Data for Q L /L, the "parallel" chiral 
glass correlation length divided by system size, as a function 
of T for different system sizes. 



which would indicate spin glass order. However, the situ- 
ation for the largest sizes is less clearcut, with the L = 32 
and 24 data only coming together at a somewhat lower 
temperature than the temperature where the L < 16 
data intersect. Furthermore the data for the largest sizes 



at the lowest temperatures does not depend strongly on 
size, indicating close to marginal behavior. This is quali- 
tatively in agreement with Campos et al.— . An enlarged 
view of the important region is shown in Fig. [6] 

Data for the parallel chiral glass correlation length is 
shown in Fig. [7] The main features are the same as found 
for the spin glass correlation length in Fig.0 The smaller 
sizes show clear intersections, but the L = 32 data lies 
lower, though perhaps not quite to the same extent as 
in Fig. [5l As for the spin glass correlation length, the 
two larger sizes L = 32 and 24 only come together at a 
lower temperature than the temperature (range) where 
the smaller sizes intersect, and do not splay out at still 
lower temperatures, at least in the range of T studied 
and within the error bars. An enlarged view is shown in 
Fig.[5J In our view, the data in Figs. [5] and [7] are not very 
different; in particular, in both figures, the data for the 
largest two sizes merge at about the same temperature 
and stick together at lower temperatures. If we look at 
the largest three sizes there is a somewhat greater ten- 
dency for the chiral data to splay out at low-T. Clearly 
larger sizes are needed to be sure of the trend in the 
thermodynamic limit. 

With this data it is difficult to come to a firm con- 
clusion about the nature of the transition. Clearly the 
smaller sizes have large corrections to scaling behavior, 
and so, at best, only the two largest sizes are in the 
asymptotic scaling region. The data is consistent with a 
line of critical points analogous to that in the KTB tran- 
sition, as proposed by Campos et al., in which the data 
for large sizes would merge at a single temperature and 
remain independent of size at lower temperatures. This is 
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a scenario in which the lower critical dimension di is equal 
to three. Another scenario consistent with our data, in 
which di is also equal to 3, is that Tsg — Tea = but 
with an exponential divergence of the correlation lengths 
at T = 0. In that case, it is likely that the data for large 
sizes would join a common curve, but the temperature at 
which the common curve is joined would decrease as the 
system size increased. This scenario is found in the two- 
dimensional Heisenberg ferromagnet, whereas the critical 
line of the KTB theory is found in the two-dimensional 
XY-ferromagnet. Our data is also consistent with the 
possibility that, for the large sizes, the curves weakly 
intersect, implying a finite Tsg & n d a lower critical di- 
mension slightly less than three. 

Unfortunately, because of the crossover effects in the 
data, it does not appear possible to give a meaningful 
estimate of the critical exponents, v and r\. 

IV. CONCLUSIONS 

In agreement with Campos et al.— and Hukushima and 
Kawamura^, we find a crossover in behavior for system 
sizes larger than about 16. Results from the larger sizes 
indicate a more marginal behavior than those from the 
smaller sizes. Unlike Kawamura and collaborators^^, 
we do not find that this effect is very different for the 
spins and chiralities. Compare, for example, Figs. [5] and 
[71 However, larger sizes would be needed to confirm that 
the asymptotic behaviors of the spin and chiral glass cor- 
relation lengths are indeed similar. Our data extends up 
to L = 32, somewhat larger than the sizes [L < 20) stud- 
ied by Hukushima and Kawamura. Our range of sizes is 
the same as that of Campos et al., but we are able to go 
down to lower temperature, in particular below the puta- 
tive spin glass transition temperature. Whereas Campos 
et al. argue in favor of a critical line below Tsg, m our 
view, other possibilities can not be ruled out, such as a 
transition at a lower value of Tsg or even an exponential 
divergence of the correlation length at T = 0. 

To distinguish between these scenarios would require a 
study of still larger sizes. It will be difficult to go to very 
much larger sizes without a better algorithm, since the 
present study used a quite substantial amount of com- 
puter time. An unfortunate feature of the present algo- 
rithm is that parallel tempering for vector spin models 
requires a large number of temperatures. The large num- 



ber arises from the temperature independent specific heat 
at low temperatures, which, in turn, comes from from 
a rather trivial feature: Gaussian spinwave fluctuations. 
However, it is difficult to see how to eliminate their effect, 
and thereby reduce the number of temperatures. 

It is interesting to note that we (and also Campos et 
al?i£) have been able to study larger sizes for the Heisen- 
berg spin glass (L = 32) than has ever been done for the 
Ising spin glass. For example, Katzgraber et ali2£ stud- 
ied the Ising spin glass using considerable CPU time, but 
were still not able to equilibrate sizes larger than L = 16 
near Tsg for the case of Gaussian interactions. With ± J 
interactions, where "multispin coding" speeds the code 
up further, they could go up to L = 24. 

That one can equilibrate larger sizes in the Heisenberg 
case is surprising bearing in mind that (i) the Heisen- 
berg algorithm is more complicated and so one sweep 
takes more CPU time than for the Ising case, and (ii) 
the transition temperature is significantly lower for the 
Heisenberg case. Noting that the mean field transition 
temperature is TgQ = s/z/m, where z is the number 
of neighbors and m the number of spin components, if 
we take Tsg = 0.145 for the Heisenberg case, the ratio 
Tsg/Tsg F is about 0.18. For the Ising case, one finds 
e.g. Ref. [H, Tsg - 0.95 so T SG /Tl G F ~ 0.39 which is 
more than twice the corresponding ratio for the Heisen- 
berg spin glass. 

The fact that, despite all this, one can study larger 
sizes in Heisenberg spin glasses than in Ising spin glasses 
indicates that barriers are smaller in the Heisenberg case. 
Another way of putting this is that the extra degrees of 
freedom in the Heisenberg model allow the system to find 
a way round barriers, which would have to be gone over 
for the Ising spin glass. 

Overall, our results indicate that spins and chiralities 
behave in a quite similar manner, and that the lower 
critical dimension of the Heisenberg spin glass is close to, 
and possibly equal to, three. 
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